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§ ■ Abstract 



We review in details a recently proposed technique to extract information 
about dynamical correlation functions of many-body hamiltonians with a few 
Lanczos iterations and without the limitation of finite size. We apply this 
technique to understand the low energy properties and the dynamical spec- 
tral weight of a simple model describing the motion of a single hole in a 



"j^ ■ quantum antiferromagnet: the t — Jz model in two spatial dimension and for 

g 

I ■ a double chain lattice. The simplicity of the model allows us a well controlled 

^ : 

^ . numerical solution, especially for the two chain case. Contrary to previous 



approximations we have found that the single hole ground state in the infinite 



^ . system is continuously connected with the Nagaoka fully polarized state for 

J2 — > 0. Analogously we have obtained an accurate determination of the dy- 
namical spectral weight relevant for photoemission experiments. For = 
an argument is given that the spectral weight vanishes at the Nagaoka energy 
faster than any power law, as supported also by a clear numerical evidence. It 
is also shown that spin charge decoupling is an exact property for a single hole 
in the Bethe lattice but does not apply to the more realistic lattices where 

the hole can describe closed loop paths. 
75.10.Jm,75.40.Mg,71.10.-Fx 
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I. INTRODUCTION 



After the discovery of the high-T^ materials there have been an increasing attention 
to the study of strongly correlated fermion systems. In his pioneer work Anderson first 
pointed out that a doped Mott insulator, described by the one band Hubbard model, should 
contain the basic physics of the high-T^ superconductors. The Hubbard model was proposed 
independently in 1963 by Gutzwiller, Hubbard and Kanamori |^-|^. Despite of its simple- 
looking hamiltonian and a lot of efforts, the physics of the model is still a subject of debate. 

In the large U limit, it is more convenient to use a canonical transformation to project 
out the doubly occupied sites, costing energy U. This leads to the so-called t-J model: 

H =-t + h.C.) + J ^ (Si • Sj - ^niTtj). (1) 

<i,j>,a <i,j> 

where the constraint of no double occupancy is understood. Here cj(cj) creates (annihilates) 
an electron at site i, rii = J2''^iu is the corresponding density operator with = cj^Ci^, 
the symbol < i,j > means summation over nearest neighbors, J = ^ is the superexchange 
coupling and finally the spin density operator Sj is defined by the Pauli matrices a: Sj = 

cr,cr' 

Starting from a more realistic three band Hubbard hamiltonian describing the Copper- 
Oxygen layer, Zhang and Rice showed that the low-energy physics of cuprate superconduc- 
tors is determined by the singlet state formed by the additional hole on oxygen with the 
existing hole copper. They also showed that the hopping of this singlet is described by 
an effective one-band t-J model, with J and \J essentially unrelated. This makes the t — J 
model even more appropriate to describe the Copper-Oxide planes at low energies rather 
than the original one-band Hubbard model. 

During the last few years, a huge amount of analytical and numerical work has been 
devoted to the study of the t — J model. |0] However the physics contained in the t-J model 
is far from clear because of the interplay between antiferromagnetic long range order and 
charge degrees of freedom. Exact numerical methods are usually limited to a very small 



2 



linear size in more than ID and practically nothing about the low energy physics has been 
understood numerically. 

In this paper we will make a detailed numerical study of a simplified version of the t-J 
model, i.e. neglecting the spin fiuctuations in the exchange term (Sj ■ Sj SfS^ in Eq. (|1|)) 
and considering only the properties of a single hole. Due to the simplification of the model, 
we are able to work directly in the infinite size system and have a satisfactory description 
for the low energy dynamics of a single hole in such a simple model of an antiferromagnet: 

H = -t (49. + h.c.) + J.Y sts^. (2) 

<i,j>,<j <id> 

A basic motivation of this work is to obtain reliable numerical results on this simple t — 
model since it may be useful for further investigation of the more interesting t — J model. 
For instance, the t — model was often used in the past to test several approximations 
on the t — J model. In fact most of the physical properties of the t-J model probably remain 
valid even for the t-Jz model. From this point of view it is worth mentioning that the t — J 
model and the t — model have the same limit for J or — and for infinite spatial 
dimensions the two model coincide, since the spin fiuctuations are irrelevant in this limit. 

m 

The single hole problem represents certainly a further simplification but is still physically 
relevant, since the single particle excitations in magnetic insulators can actually be studied 
by the photoemission and the inverse photoemission spectroscopy. 

For a single hole a rigorous theorem proven by Nagaoka is known for = 0. The so 
called Nagaoka theorem states that for a bipartite finite lattice in more than one dimension 
{d > 1) the ferromagnetic state with maximum spin is the unique ground state with energy 
Cat = —zt., in any subspace with given total spin projection. 

A complete description of the one hole spectrum in the = limit, not limited to the 



ground state, was first given in where the so-called "retraceable path approximation" 
(RPA) was introduced. In the Ising limit as a hole hops in a Neel state, it scrambles the spins 
along its path. In order to return the spin configuration to its original state, it was argued 



that, to a good approximation, one can consider only paths in which the hole retraces its 
path back to the origin, thereby returning all of the spins to their original position. In this 
approximation one can write down an explicit analytic solution for the Green's function: 



G^i-) = ^ 1. J (3) 



and the spectral weight A{k,u) = —^^^IrnGkiyj) reads 



A{uj) = \ \ „ . (4) 

where z = 2d is the number of nearest neighbors. The spectral weight is completely inco- 
herent and dispersionless with a one particle band, which is 75% narrower (in 3D) than the 
noninteracting band. The RPA is exact in ID and for the Bethe lattice (where Nagaoka 
theorem does not apply), and recently it has been shown to be exact in the limit of infinite 
spatial dimensionality. |T^ 

For finite J^/t, the competition between the kinetic energy t favoring the ferromagnetic 
configuration and the exchange energy Jz favoring the antiferromagnetic alignment of the 
neighboring spins makes the problem of particular interest. By neglecting closed loops. 



e.g. in the Bethe lattice case, it is possible to derive a closed solution [|12| , p]3i , which we 
will refer in the following as the "string picture". In the string picture, the hole moves in 
an antiferromagnetic spin background, leaving behind a string of overturned spins, which 
costs an energy proportional to the length of the path. The overturned spins behave like an 
effective linear potential for the hole. In the continuum limit, valid for —>■ 0, the problem 
is reduced to a one dimensional Schrodinger equation with a linear potential V ~ xJ^ where 
X is the length of the string. 

. JAz-2) 



H = -^T^t— + '-X - 2V7^t. (5) 

ax'^ 2 

The solution of this hamiltonian leads to a series of bound states. The hole is essentially 
localized with a discrete energy spectrum: 

E„ = ar^VT^ti i'^' ' '^} f'^ - 2V^t (6) 

Zy Z — Lt 
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where a„ are the zeros of the Airy function Ai{z). Recently, it has been shown that the 
string picture is exact up to order ^, where d is the spatial dimensionality. 



Contrary to the = case, at finite the spectral weight has 5— function peaks at 
energies En. The weight of the 6— function at the lowest energy is called the quasiparticle 
weight Z and is found to vanish linearly in J^, using the continuous limit (^. |^ 

Later Kane, Lee and Read |^ introduced a self-consistent Born approximation that can 
be considered an extension of the RPA to the more physical t — J model. Their approach 
is widely accepted since they were able to reproduce the asymptotic behavior of the t — 
model in the string picture and a large amount of numerical work on small lattices seems to 
be in qualitative agreement with the predictions of this theory. |ll^-|l6| 



In the retraceable path approximation or in the string picture and similarly within the 
self-consistent Born approximation, closed loop paths are neglected in order to simplify the 
analytic structure of the one hole Green's function. However, this certainly introduce an 
approximation even in the simple t — J^ model. In fact the most important effect due to the 
inclusion of closed loop paths in this model was first noted by Trugman. In 2D a hole can 
hop around a square plaquette one and half times without disturbing the spin background 
with a net translation to the next-nearest-neighbor along the diagonal. This means that 
the hole can "unwind" the string and self-generate a next-nearest-neighbor hopping. The 
full localization of charge carriers found within the retraceable path approximation and the 
string picture is an artifact of the approximation. 

Another interesting issue is whether the spin and charge degrees of freedom are decoupled 
or not. [p!7| , p!8| Spin charge decoupling is a well known phenomenon in one dimension, where 



a one electron excitation can be decomposed into a spinon excitation which carries spin 
but no charge and a holon excitation which carries charge but no spin. We will show in 
the present paper that spin charge decoupling is an exact property of the t-J^ model in the 
Bethe lattice with arbitrary coordination number z at single hole doping and present some 
numerical work, ruling out the possibility of spin charge decoupling in physically relevant 
lattices where closed loops are allowed. 



In this paper, we make use of the Lanczos method to calculate the single hole spectral 
function in the infinite system. We will first introduce a transformation to eliminate the 
charge degree of freedom to get an effective spin Hamiltonian in a given momentum subspace 
(Section 2). Then we will use the Lanczos method in order to diagonalize the effective spin 
Hamiltonian in the infinite lattice (Section 3). We will concentrate on two chains (2C) 
and two-dimensional (2D) square lattices. The 2C case is much easier for numerical study 
compared to the 2D lattice. The 2C lattice is not a trivial lattice for the t-J^ model, because 
it has the basic properties of an higher dimensional lattice, such as the existence of closed 
loop paths and the validity of the Nagaoka theorem. 

Most of the above results were obtained by use of the Lanczos Spectra Decoding method. 
This method allows us to analyze the Lanczos data in the infinite system, where only a small 
number of Lanczos iterations is possible (Section 3). The validity of the Lanczos Spectra 
Decoding is proved analytically and numerically on the Bethe lattice (Section 4). In section 
5, we will show our numerical results both for the = and finite case. Finally we will 
discuss about the formation of a ferromagnetic polaron around the hole for small J^, and 
about the analytic form of the spectral weight close to the Nagaoka energy. 



II. FORMALISM 

A. Effective Spin Hamiltonian 

The t — Jz Hamiltonian (|^) is translational invariant, so the most general one- hole state 
with total spin Sj = Mz + | and total momentum — p ( hole momentum p ) can be written: 

IV'P > = ^ E e'^^Cfi,.T«|^o > (7) 

where \So > is a pure spin state that satisfies: i) nR\So >= \So > for all sites R, ii) 
'^j.ol'S'o > = \So > i-e. the spin at the origin is fixed to (Tq = — |, iii) the total spin along 
the z-axis is a well defined quantum number in the state \So >, i-e. S'^\So >= Mz\So >■ 
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The operator in (|^) is the translation operator that brings the origin O to the lattice 
point R. It is formally defined by: 

TrCr'^^T-R = Cr+r/^ct (8) 

States with definite momentum {ipp > are eigenstates of Tr such that Tr\iPp > = e~'^^^\ipp >. 
Note also that the sum over a in the definition of {tpp > in (0) is used only for later conve- 
nience. 

A pure spin Hamiltonian can be derived as following. We evaluate the expectation value 
of the t-Jz Hamiltonian on the one-hole state \tpp >, Ep = < ipplHlipp >, and consider first 
the case = 0: 

= ^ E e-P(«-«^) < So\T^rA,,^P CR^r,,A,.CR,,.,TR,\So > (9) 

-Rl ,R2,R,Tf_i ,a,ai ,(T2 

Now i) (Ti = (72 since the total has to be conserved in order to give a non-vanishing 
contribution in Ep. ii) R = R2 and R + = Ri, otherwise we create a doubly occupied site 
which is projected out either by P or by overlapping with < So\- Finally we obtain that 
Ep =<So\H^^^So > with: 

Hf=t E e-'^-^T^rM.yCr.,<^cl^co,.'). (10) 

ifp-^-^ acts only on spin states. In fact for S = | it is possible to express the term between 
parenthesis in terms of spin operators only: 

XR^,Ri = E ck,'7'^R^,<^CR„<7CR„a' = ^riRMR^ + 2S R^ " S R^ (ll) 

Finally Hp^'l' reads: 

H;" = tY.e-''^'T.^.Xo,r, = tT.^''^'Xo,rT^,- (12) 

Note that the product of the operators T_t-^ and Xo,t^ appearing in -ffp-^-^ permutes spins 
and leaves the origin O unchanged. Moreover, since the operator x interchanges the spins 
at sites Ri and Rj, we have that XRi,Rj = XRj,Rr and that X%,Rj = 1- 
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Analogoously we can extend to J2 7^ the derivation of a pure spin Hamiltonian and 
obtain: 



H;" = tY.e'''^Xo,rJr, + J.i: S%^S%^ (13) 

The full Hamiltonian commutes with So (since it actually does not depend on the spin at 
the origin) and thus it can be defined on — 1 sites. 

In conclusion for any eigenstate 15*0 > of H^^^ with definite Sq\So >= —^\So > we 
have an eigenstate of the t — Hamiltonian written in the form (^. In fact, by use of 
the variational principle any eigenstate of (§) or (|1^) is obtained by ^ < ipp\H\ipp >= 0, 
< So\H\So >= with the condition < -ippl^p >= 1, < So\So >= 1 respectively. Since 
Ep =< ipp\H\ipp >=< So\Hp^^So > and < ipplipp >=< So\So >, it clearly follows that all 
the eigenstates of Hp^^ with definite spin at the origin define a true eigenstate of H by use 
of (|^. In conclusion the one hole problem is mapped onto a diagonalization of a pure spin 
Hamiltonian H;ff for the given momentum p. [|l^] A similar effective Hamiltonian can be 



obtained for the t-J model, by substituting 5*^ with S in (|T3D. 

It is interesting to note that Hp-^-l" is not translational invariant and that the momentum 
of the hole appears as a simple parameter. It is just this property that allow us to diagonalize 
the t-Jz model in certain momentum subspace in the infinite lattice. 

As it is shown in the Appendix many useful dynamical quantities such as the Green's 
function and current operators can be easily translated in terms of spin operators acting on 
this spin space in which the hole is fixed at the origin. 

B. Green's Function for X = 



Following Brinkman and Rice ||Tl|], we can expand the Green's function G in terms of the 
momenta < H\{HpffY\H > of the Hamiltonian on the translation invariant ground state of 
the undoped system \H >. For = 0, at vanishing doping, the hamiltonian is classical and 
I if > is given by: 
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\H >= ^{\N > +\N' >) (14) 

where \N > and |A^' > are the two possible determination of the Neel state. 

The one hole Green's function can be generally written as a summation of all possible 
paths traced by the hole during its motion on the lattice. A path is then defined by a set of 
coordinates {-R/}" with / = 0, n, which are connected by nearest neighbor vectors 

Ri-Ri-i = T^,. (15) 

Among all the possible paths it is useful to identify the skeleton ones. ||2^ A path {Ri}"' of 
length n is a skeleton path if -R/+i 7^ Ri~i for any I = 1, ...,n — 1. By definition, for any 
skeleton path, the hole never retraces its path immediately. It is clear that all the remaining 
paths can be obtained by dressing each site Ri of the skeleton path by all possible retraceable 
paths. The retraceable paths can be then summed exactly. 

Instead of giving the detailed derivation, we will present the final result and discuss some 
of the special cases. The most general Green's function can be written in a formal expansion 
of skeleton paths of length 2n: 

G{R,u) = GBR{uj)Y,K{ujf^C^{R) (16) 

n 

where Gbr is the Brinkman Rice result (^, which includes only retraceable path contribution 
starting from the origin R = O and coming back to the same site: 



while the function K{uj) is the exact contribution of all possible retraceable paths on each 
site of the skeleton path different from the origin: 
t 1 1 



cjl-S^(cj) 2{z-l)t 



uj — sgn(u;)i/a;2 — 4(z — 1)^^ 



(18) 



I 2n^ 



Finally in (|T^ the coefficients C„(-R) are determined by: C„(-R) = ^{{Ri}' 

all skeleton {i?;}^" 

where !]({/?/}") are spin correlation functions, which, generally speaking, depend on the 
path. 



^{{RlV) = < H\xR„,R„.,XRr.-uRu-2 ■ ■ ■ XRr,Ro\H > (19) 

In the Ising case Q{{Ri}"') vanish for odd n and are simply equal to either zero or one 
for even n depending on whether the skeleton path change the Neel order or not. Instead 
they are complicated correlation functions in the t — J model for J ^ 0. In this case the 
expression (|TB|) is still valid provided we include the odd n contributions. 

Two possible limits can be exactly solved using the previous expression (|16|) for the 
Green's function. 

1. Bethe Lattice Case 

The Bethe lattice is defined on a Cailey-tree with coordination z and for z = 2 it coincides 
with the one dimensional chain. In the Bethe lattice there is only one path that connects 
two arbitrary sites of the lattice because there are no closed loops in this lattice (see e.g. 
||T6| for a more detailed definition of the Bethe lattice). As in one dimension there is only 
one skeleton path connecting the origin to a given point R and Q{{Ri}'^) = Q{Rn), which 
is only a function of the final position Rn. We can immediately write the exact expression 
of the Green's function on a Bethe lattice using the previous general expression (|T6|): 

G{Rr„u) = n{Rn) G^'^'^^iu) (20) 

where G^^'^'^ = Gbr{uj) K{uj)"'. It is interesting that the "strong correlation" in the one hole 
problem in the Bethe lattice is only contained in the static function Q{Rn), since G^'^^^{uj) is 
exactly the free electron Green's function. This has direct consequences on the spin charge 
decoupling. (See next section) 

2. t-J^: Brinkman-Rice Almost "Exact" 

Let us consider the diagonal Green's function G{R = 0,uj) for the t-Jz model. In 
the general expression (|l^), Q{{Ri}"') is either one or zero depending on whether a given 
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permutation preserves the Neel order or not. The shortest skeleton paths with f2({i?;}") = 1 
is known to be three times the path around the elementary square plaquette. Thus n = 12 
for such a skeleton path. There are 8 possibilities to build such a path starting from the 
origin in 2D ( 4 neighbors times 2 possible opposite directions) and 2d{2d — 2) in dimension 
d. Then the next leading correction to the Brinkman and Rice result is 

G{R = 0,uj) = Gbr{oj) (l + 2d{2d - 2) K{u)^'^ + ...) (21) 

Since |-ft'(co')| < ^^d-i ^^^^ ^' correction to the Brinkman and Rice result turns out 
to be less then 1.1% in d = 2 and 0.15% in c? = 3. This has also been noted in ^ expansion 
T0| but appears also quite natural in this formalism. 



C. Spin Charge Decoupling 

In ID, if spin charge decoupling occurs, the one particle Green's function can be written 
as a simple product of a spinon contribution and a holon contribution: 

G{R, t) = t) G''°^""iR, t) (22) 

This property is exact in = 1 for the infinite U Hubbard model and is asymptotically valid 
for large (i?, t) at finite U . 



The spin charge decoupling manifests itself in the one particle Green's function and in 
principle can be detected even at higher dimensionahty as speculated by several authors 
following P.W. Anderson. 0] For J, — 0, or f/ ^ oo one expects no dynamics for the 
spinous and that the holon contribution has exactly the free particle behavior, because there 
is only a unitary charge carried by the single hole. 

Thus, as a consequence of spin charge separation, the one hole Green's function should 
be written in the following way: 

G{R,t) = n{R) G^''\R,t) (23) 

The free electron Green's function is nothing but the free propagator in the Nagaoka limit: 
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G^'"'(k,u) = (24) 

uj — tk + tosgnu 



where = 2t{coskx + cos ky) is the energy of a free hole. 



By Fourier transforming in time and by taking the imaginary part of the equation (23) 



we obtain the spectral weight as a function of the final position of the hole and the frequency 



A{u, R) = Q{R) A'''\u, R) (25) 

The previous expression is exact in one dimension even for this simplified t — Jz model 
where the spinon function is particularly simple Q{R) = 5r^o- [lU By Eq. (|20|) spin-charge 
decoupling is valid in the Bethe lattice for the single hole problem. 

Equation (^) can be considered as a direct and measurable consequence of spin-charge 
separation. In fact , by measuring the spectral weight for two different positions of the hole 
we should get that the ratio: 

A{uo,R') n{R') 



A{lj,R) Q{R) 



independent of uj (26) 



The ratio 4^^f^, according to expression (PBD, should be independent of uj if spin-charge 
decoupling occurs. This is unlikely in dimension higher than one as the presence of the 
skeleton paths strongly renormalizes the spectral weight with the distance of the hole from 
the origin. p2[ The absence of spin-charge decoupling is particularly evident already for the 
2C case (Fig.p, where A{u, R)/A{uj, R = O) is u dependent for relatively large R and quite 
low energy. 

Based on the numerical and the analytical results, we conclude that spin-charge decou- 
pling can occur mainly in lattices where closed loop paths are forbidden by the geometry. 
In these type of lattices the Nagaoka theorem cannot be applied. We have thus found an 
interesting relation between the Nagaoka theorem, spin charge decoupling and presence or 
absence of skeleton paths in a given lattice. Of course we cannot rule out more complicated 
form of spin-charge decoupling such as the one discussed in [^] , or that at finite doping the 
situation may change, fl^ 
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III. LANCZOS SCHEME 



The Lanczos technique is widely used in strongly correlated electron systems. Contrary 
to the quantum Monte Carlo technique, it does not suffer the "fermion sign problem" or any 
other instabilities at low temperature. Dynamical correlations can be easily obtained using 
this technique. However, it has been restricted so far to small systems, typically 4x4 (for 
Hubbard models), or at most 26 sites (for t-J and t-J^ models) and 36 sites (for Heisenberg 
models). [T^ On the other hand a systematic way for a finite-size scaling analysis in doped 
system such as the t-Jz model is not known yet and some infinite-system properties are still 
unclear or even misleading. 

In this section, we will develop a scheme, which allows us to analyze the infinite system 
Lanczos data in a very efficient way, so that we are able to calculate the spectral function 
of the t-Jz model with good accuracy. 

The Lanczos method is devised to diagonalize huge Hamiltonian matrix with dimen- 
sion Nh- The method starts with a trial wave function 0^. A new basis is generated by 
Hamiltonian multiplication, 

Si = W\(t)T> for z = 0,l,...,n (27) 

Then an orthogonal basis { } can be iteratively calculated, after orthogonalization of the 
vectors Sj. Formally we have 

6j+i|ei+i> = H\ei > -ai\e^ > -&j|ei_i > 
ai = < ei\H\ei > 
bi+i = < e^+i\H\ei> (28) 

where Bq = 0, |eo >= \<f)T > ■ In the Lanczos basis, the Hamiltonian turns to a tridiagonal 
form where are the diagonal elements and bi the off-diagonal ones for i = 1, . . .n. 

For n = Nh, the spectrum of the tridiagonal matrix coincides with the one of the original 
Hamiltonian. Unfortunately the dimension of the Hilbert space is usually given by ~ 
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and it is prohibitive to perform a full diagonalization with the available computers. 
Nevertheless the ground state and the corresponding energy converge for relatively small 
n ~ 10^ << Nh- This justifies the success of the method which enables to restrict the 
diagonalization to a very small basis n « Nh where, according to the Ritz theorem, the 
variational principle applies for all the eigenvalues of the smaller tridiagonal matrix and in 
particular for the ground state energy. 

A. Lanczos scheme in an Infinite Lattice 

In an infinite system, the Lanczos scheme (^81) can be applied efficiently, provided there 
is a proper way to define a simple finite basis to represent the vectors s„ = H"'\tpT > (p7|), 
generated by the iterative application of H to the trial state. This is in fact the case for 
the effective t-J^ Hamiltonian (|13|), where Sj is represented by overturned spins on the Neel 
state localized around the hole. 

For any fixed momentum p, we start from a Neel state with a hole at origin. Since 
the Jz term of the Hamiltonian is diagonal, the only part of the effective t-J^ Hamiltonian, 
relevant to generate new states, is the kinetic term. In each multiplication of the effective 
Hamiltonian, the hole is translated to its z nearest neighbors by the translation operator 
Tt-^ (see Fig.(|2D). Then the spin exchange operators Xo,t^ move the hole back to the origin, 
leaving an overturned spin background and generating z new states. 

The possibility to work with a finite basis even in the infinite system was first noted by 
Trugman [Q. In fact the overturned spins are located within a region around the hole with 
radius n. We can thus update only the defects over the Neel state, which are finite at any 
finite number of multiplications of H. After n steps, the Hilbert space is finite having at 
most dimension z^. 

This exponential growth of the Hilbert space ~ makes the problem intractable even 
for relatively small n. Fortunately many of the generated states appear several times during 
the expansion of the Hilbert space, due to the presence of the Trugman-like paths, |g] and 
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also due to the translation symmetry implicitly exploited by use of the effective Hamiltonian 



(p!3|). After all, the dimension of the Hilbert space turns out to be considerably smaller than 
the previous estimate, and in fact it grows much slower than 2;", e.g. ~ 1.9" for z = 3. In 
this way we have reached n = 2Q for the 2C case and n = lA for the 2D case with an Hilbert 
space dimension at most equal to ~ 12.2 x 10^. 

The smallest finite lattice which contains the full information of the first exact n Lanczos 
steps has linear dimension 2n + 1. Therefore, our results correspond to a 53 x 2 lattice in 
the 2C case and to a 27 x 27 in 2D, which is by far larger than the size of a typical finite size 



Lanczos calculation [Q. However the finiteness of n often leads to appreciable systematic 
errors. In these cases we overcome this difficulty by carrying out a systematic extrapolation 
in l/n — s> 0, following a scheme which is analogous to the finite size scaling analysis for finite 
lattice calculations. 

We conclude this section with some technical comments about the algorithm. The basis 
generated by the iterative application of H^'^-l' does not depend on the momentum p of the 
hole. Hence we only need to generate it once, which takes less than two hundred seconds 
of CPU time on a Cray-C90. After that, we can do the usual Lanczos iterations for fixed 
parameters p and J^, which typically takes 10^ -r- 10^ seconds of CPU time on the same 
computer. 

B. Lanczos Spectra Decoding 

With the n + 1 eigenvalues Ei and eigenfunction > of the Lanczos matrix truncated 
after n step, the spectral weight can be formally calculated as 

Aik,u;) = Imi(vl>^|^^|vl/^) 

(29) 

In the following we assume that the energies Ei are set in ascending order : -Ei+i > E^. 

As a result of the finiteness of the restricted Hilbert space, we get a sum of (5-functions 
in the spectral weight at any fixed n. Due to the Lehmann representation of the spectral 
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function this feature is also present in any finite size calculation. In this case the 



thermodynamic limit is obtained by smoothing the 5-functions in Eq. with Lorenzians 



of a given small width 5, [TIJT 



5{u;-E,) ^ Im (30) 

uj — Li — to 

and then taking the limit 5 — > 0. For small finite S, reasonable results can be obtained by 
the finite size Lanczos algorithm, provided that the resolution of the energy levels becomes 



much smaller than S, i.e. n large enough but still much less than Nh- |114[ 

Numerically we have been able to perform n = 26 and n = 14 Lanczos iterations for 
the 2C and the 2D lattice respectively. Even though the ground state energy is already 
well convergent, such a small number of Lanczos steps is usually far from enough for a 
good estimate of the spectral weight. In fact by the conventional method of smoothing the 
5— functions described in (pOD, either one miss the details of the spectral weight for large S 
or in the opposite case one gets too rapid oscillations which are obviously unphysical. 

A more efficient method for evaluating the spectral weight was recently introduced by us. 
In the following, for reason of completeness we will give a brief review of this new method 



named Lanczos Spectra Decoding. In this simple method, we introduced an interpretation 
of the Lanczos scheme. With this interpretation, the spectral function can be calculated 
accurately, efficiently and easily even with a small number of iterations n. As we have 
seen in fact the Lanczos scheme in the infinite system has a computational cost growing 
exponentially with n, whereas in any finite size calculation the algorithm is only linear in 
n. We expect therefore that our method is essential in the first case but maybe helpful only 
for a small computer-time factor in a finite size calculation. 

As well known the spectral weight A{uj) is a distribution that may be divided into two 
parts A{uj) = Acoh{^) + ^mco/il'^), a coherent one Acoh{^) which contains only 5— function 
contributions and an incoherent one Amcohi^^) which is a continuous and usually smooth 
function of u. The Lanczos spectra decoding exploits the smoothness properties of Aincohi^^) 
in a simple and efficient way. In the following we therefore assume that the spectral weight 
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is incoherent. This is not a hmitation since coherent peaks can be easily separated out from 
the incoherent part by identifying all the quasiparticle weights — |(\E'j|\E'r)p that remain 
finite for n — > oo. 

n 

If the spectral weight is incoherent , since by completeness ^ Zi — \^ one expects that 

i=0 

(X ^. Thus 

Z{u) = (n + \)Zi (31) 

may define a smooth function of lo at the discrete Lanczos energies cu — Ei. The full spectral 
weight is then closely related to the previous function: 

A{u) = Z{u)pl{ou) (32) 

where pl represents the Lanczos density of states (LDOS) in the restricted Hilbert space 
generated by the Lanczos algorithm: 

1 

Pl{u:) = ^Y.^{u - E,) (33) 

oo 

where the factor is determined by the normalization condition / depL{e) = 1. For 

— oo 

n — Nfi in a finite system, p^ coincides with the actual density of states of the many body 
system, but, in an infinite lattice { — oo) , this is not generally true, as it will be shown 
for the Bethe lattice case. 

By definition, the number of states dN between energies e and e + de is given by 

dN ={n + l)pL{e)de. (34) 

So the Lanczos density of states can be calculated as 

I dN 

Finally , using finite difference instead of differential, the coarse grained Lanczos density of 
states can be estimated up to order O(^) by 

''^'^ = '''' 
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where the energies 



he at the middle of two consecutive eigenvalues. 

The function Z{uj) which is known at energies Ei can be easily interpolated at the energies 
Q where also pz, is known: 



If Z{u!) is a twice-differentiable function, eq. ( P^ ) is accurate up to 0{^) as well. Thus 
within the same accuracy A{uj) = Z{lj)pl{uj) easily follows at the discrete energies u = ef. 



The knowledge of A{u) close to the previous "special" points q can be easily achieved 
by standard piecewise interpolations or extrapolations. At the end, we can verify the sum 
rule / A{uj)duj = 1, as a check for the accuracy of the calculation. 

IV. LANCZOS SPECTRA DECODING ON THE BETHE LATTICE 

In this section we will show that Z{uj) and Pl{'^) are well defined functions and can be 
calculated exactly in the Bethe lattice. 

On the Bethe lattice, the problem is exactly solvable because the skeleton paths are 
absent and the retraceable paths can be summed analytically. This exact solution has two 
meaning for us: (1) as a test for our scheme and assumptions, (2) as an hint to interpret 
the Lanczos scheme for the physical 2C and 2D lattices. 

Using the Lanczos basis defined in Eq. (pHf) , the Hamiltonian has vanishing diagonal part 
(a„ = 0) and 



Z{t,) = + Zi)/2. 



(38) 




(39) 




(40) 



n 



for n > 1 
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This easily follows since, on the Bethe lattice, each multiplication of the Hamiltonian gen- 
erates z — 1 new states, except for the first iteration, which generates z states. Due to the 
simplicity of such Lanczos matrix it is then possible to compute analytically Z and pL in 
eqs.(0,g|). 

In fact the wavefunction components i = 0, 1, . . . , n of an eigenstate with energy uj 
satisfy the iterative relation in the Lanczos basis: 

y/zipi = ^i^o i = 
\^ipo + Vz -I1P2 = ooipi i = 1 
{ipi+i + iJi-i)y/z - 1 = uipi i>2 

From the third equation solutions are possible for lo"^ < t^j^ where eBR = -2ty/I^ is 
the Brinkman-Rice approximation for the single hole ground state energy. Then the com- 
ponents of the eigenstates on the Lanczos basis are given by: ipi = lie (^A ~ , where 

Jlu^ — — UJ 

A = — and the complex number A is obtained by the first two equations, 

oj"^ — z 

yielding ipi = ^ "ipo and 1^2 = , = ipo- Finally ipo is determined by the normalization 

^ ^ziz - 1) 

condition Yl'^f = I5 and the function Z{uj) = lim ipQ{n + 1) reads: 



The Lanczos density of states (|33D is then determined by: 

Pl{u) = A{u)/Z{u) (42) 

where A{lj) is the Bethe lattice density of states in Eq.(^). In this case, Pl{^jj) has singular 
inverse square root behavior near the band tail. 

By use of the Lanczos Spectra decoding introduced in the previous section we have 
obtained a good approximation to the exact solution with n ~ 10, much less than in the 
conventional calculation. ||2H 



At finite J^, the only change is with the diagonal part of the Lanczos matrix 
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ao = -^z (43) 
= :^(32-l + 2(^-2)(z-l)) for 2^0 (44) 

This result can be easily obtained by counting the number of broken bonds in the states 
generated at different Lanczos iterations. In this formalism we obtain therefore that the 
motion of a hole in a Bethe lattice is exactly equivalent to a one dimensional motion of a 
particle in a linear potential, provided we identify the distance of the particle from the origin 
with the label i of the Lanczos basis. 

By diagonalizing numerically the Lanczos matrix for large n we can easily confirm the 
prediction of the long wavelength hamiltonian (|^), namely that the spectral function is 
fc— independent, and that A(uj) contains only 5— function peaks. In fact due to the linear 
potential all of the one-hole states are localized states. It is important to remark that 
the corrections to the asymptotic — behavior is quite important for the quasiparticle 
weight Z even for very small Jz- For Z the correct linear behavior oc Jz is found only for 
Jz ~ 10^^ (see Fig.(^)), whereas the behavior of the ground state energy and the gap are 
more reasonable. 



V. RESULTS AND DISCUSSION 



A. Spectral weight for the Jz = case 



For = the spectral weight A{uj) is an even function of uo and in the following we 
will concentrate on its negative frequency region. In this model we found p5[ a sharp 



peak in the spectral weight located at an energy close to the retraceable path prediction, 
esR = —2t\/z — 1, and, a second peak at energy ~ —t. In the 2D case, the spectral weight 
looks similar, although the first peak is rather small. In ID, the exact BR solution leads 
only to one peak but with a divergent spectral weight ~ ^J~eB bottom of the 

band. Already in the 2C case such a divergence disappears within the retraceable path 
approximation, as well as in our numerical scheme, which includes all closed-loop paths. 

20 



We also found that the first peak in the spectral function has a remarkable dispersive 
feature although the bottom of the spectrum appears fc-independent. The dispersion of the 
first peak is not present neither in ID or in infinite dimension and the importance to 
go beyond the retraceable path approximation is already clear even in 2D. 

For the density of states D{uj) = J -^yA{uj,p) = A{R = 0,uj), our results present 
some small oscillations around the retraceable path analytic solution @. In 2D however the 
retraceable path expression for z = 4 seems already quite accurate, at least away from the 
band tails. ||25| All the above results have been confirmed recently . |^ 



As discussed in the introduction a key question is the determination of the band edge 
energy e^, i.e. the threshold energy where the spectral weight begin to vanish. As a first 
step we identify the Lanczos ground state energy for n — oo as e^, which may or may 
not converge to the Nagaoka energy ejy. For instance, by neglecting closed loop paths, one 
obtains the Brinkman-Rice energy es/j as a variational estimate of Eoo- In the retraceable 
path approximation coincides with Eoo, as one generally expects. 

In order to have an accurate estimate of E^^, it is useful to have a guess about the 
asymptotic behavior of the quantity A„ = En — E^o for n oo. The way A„ vanishes 
for n ^ oo is related to the form of the Lanczos density of states at low energy. In the 
Brinkman-Rice case the exact solution in (§) gives PL(e) ~ — cb)"^''^- Thus, using 
Eqs. (|36D , A^^/^ ~ yielding A„ ~ For = the inclusion of the skeleton paths 
seems to support a finite Lanczos density of states (see Fig.(^)), yielding, by the same 
argument A„ ~ ^. We have plotted in Fig.(||) the estimated ground state energies as a 
function of 1/n for several momenta for the 2D and the 2C cases. Many of the estimated 
Lanczos energies - exact upper bound of the true ground state energy- are clearly below esR 
(even for the 2D case also shown in the picture). Thus, a previous suggestion that the one 
hole energy in a quantum antiferromagnet should be close to eBR p6|p7| is not confirmed 



by our numerical results. In Fig.(|^) -2C case- it is remarkable that all of the extrapolated 
energies are very close to the Nagaoka energy, independent of the momentum of the hole, 
i.e. = -3 ±0.02. 
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The above results give a robust evidence that Pldos is finite up to the Nagaoka energy. 
Even in this case the spectral weight A{uj) = Z{uj)pL£,os{i^) can in principle vanish for 
Lu > eN due to the vanishing of the factor Z{uj). This is the scenario suggested in |^ using 
the expansion for the Green's function in terms of K{uj) (p!6|). In the interval K{uj)a < 1 
(for \uj\/t > -^^^ + « ) when the expansion (^) converges the Green's function is surely 



real, thus determining a lower bound for e^. After some extrapolation it was found in ||2^ 



that a < z — 1, i.e. should be higher than for the 2C or the 2D case using the first 
18 or 12 coefficients Cn{R) of the Green's function expansion, respectively. Although the 
above analysis is surely correct, the basic conclusion is affected by systematic errors due to 
the limited knowledge of a few coefficients in the expansion. 

As it is shown in Fig.^, this scenario looks unlikely within the Lanczos Spectra decoding 
method because Z{ijj) seems to be smoothly connected to the Nagaoka energy. It is clear 
however that this is not enough for a definite conclusion. 

In order to solve the latter controversy without relying on the Lanczos spectra decoding, 
we have reproduced the Miiller-Hartmann- Ventura expansion (Table I) and extended it up 
to the first 26 coefficients for the 2C case. Since we have generated all of the Sj, defined in 
(P7|), up to i = n = 26 for the 2C case, the total number of paths of length 2n can be easily 
calculated as 

M„ =< s„|s„ > (45) 

By proper undressing the retraceable paths, we can convert the total number of paths M„ 
to the number of the skeleton paths C„. 

As it is shown in table I the Miiller-Hartmann- Ventura extrapolation: 

Cn{R) = CiR)j^ (46) 

used to determine the radius of convergence a is not stable when large skeleton paths are 
included. For large n, we find that a, (3{R = O) and C{R = O) are always going up. At 
n = 26, C{R = O) becomes 10 times larger than their value obtained for n = 18, while a 
changes from 1.88 to 1.91. (Table I) 
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Instead of using the fit (|46|), we apply tlie well established ratio method, which is well 
known in the study of critical phenomena. Using this method we evaluate the radius of 
convergence a, and the power law exponent 9 describing the vanishing of the Green's function 
at the band tails, corresponding to the critical temperature and to the conventional critical 
exponent in the language of critical phenomena, respectively. [^] We define 



- (47) 



and the linear intercept at - — with two next consecutive point 



s 



fi{n,n- 2) = J[n/i(n) - (n - 2)/x(n - 2)]. (48) 



In fact \i is expected to behave as 



ix{n) = a[l + ^ + 0(1)] (49) 



n 



where a gives the position of the singularity and 6 = —{1+g) gives the exponent of the 
Green's function at the band tail G{uj) oc (c<j — e^)^. This kind of analysis is much more stable 
than the extrapolation used in |2^. In Fig.(|^ the data for the linear intercept /i(n, n — 2) 



suggests that the ^ corrections are well behaved, and approximately linear for n > 20. By 
accident this number is very close to the maximum n available in the previous analysis pO| , 



giving further evidence that for this problem a very large number of coefficients C„ are 
necessary to obtain reasonably converged results. 

Then, by a linear fit in the region n > 20, we obtain a = z — 1 within one per cent, 
both for i? = and for p = 0, as an independent check that a should not depend on R 
{Cn{p = 0) = Cn{R))- This is a clear evidence that the band edge energy coincides with 
the Nagaoka energy. 

For the 2D case the data shown in Fig. (§) indicate that the above analysis is poorly 
converged due to a too much small value of n = 14. This also explains why the extrapolated 
n — i> oo results for E^o, shown in Fig(B|), does not converge to eA? in 2D. 



23 



For the critical exponent 6 we have a much less accurate result probably because the 

vanishing of the spectral weight close to ejv is characterized by overall essential singularities 

1 

of the type e """^v , as we have suggested in our previous paper. ||25[ 



Here we can give a very simple argument supporting the expected singular behavior for 
Z{u). For the spectral weight the behavior close to should be the same, because , as 
we have previously shown, the Lanczos density of states is approximately constant in this 
region. 

After n Lanczos steps the hole forms a polaron state \P > with maximum spin but total 



Sz = [^] in a region of linear size ^ ~ ra, as it is also confirmed by direct calculation of 
the spin arrangement around the hole (see next section). The overlap of the Neel state and 
this trial state is then approximately given by p9|l: 



V/2) 

where V is the volume of the polaron equal to 2^ for the 2C case and oc ^'^ for general spatial 
dimension d > 1. The quantity Zq x n according to the Lanczos spectra decoding method 
characterizes the behavior of the smooth function Z{u!) at an energy u = + const. /n. 
Solving for n from the latter equation and assuming, as we have already mentioned, that 
^ oc n we can substitute oc (u; — e^)^'^ in (|50| ) and obtain: 

A{uj) oc Z{uj) (x{uj- ejv)"^/^e"(^^ for 2C (51) 
A{uj)ocZ{uj)(x{uj-eN)"^^'^~^e~^^^ for ci > 1 (52) 

where A is an overall constant depending on the dimensionality. Using the above formulas 
we have obtained good agreement with numerical results over a range for Z covering up to 
two decades (see Fig.(||)). 

B. Spectral weight for finite 



For finite Jz a coherent part shows up in the spectral weight. However, contrary to the 
string picture, only the first few energy levels contribute to the spectral weight with true 
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S— functions (See Fig.([TO|)). In order to identify these S— function contributions we can 
check whether the quasiparticle weight Zi converge to some finite value for n oo, whereas 
if the energy level Ei belongs to the incoherent part Z{u!) = {n + l)Zi remains finite for 
n ^ oo. Another method to distinguish the coherent part from the incoherent one is to 
analyze directly the wavefunction components of the eigenstate \Efj on the Lanczos basis {ej}. 
As we have seen in the previous section the label j measures the length of the overturned 
spins in the state i. The quasiparticle weights Zi are finite only if the Neel state ej=o will 
have a non vanishing component with the state \E'j. This obviously occurs if the eigenstate 
"^i is "localized" in the Lanczos basis even for n — > oo, otherwise only a probability ~ to 
be in the Neel state is expected. 

Using the above criteria, shown in Fig.(|TT]), we have a clear evidence of a single quasi- 
particle weight for not too large, and an incoherent part which is rather similar to the 
Jz = one. For larger the incoherent part moves quite fast to higher energies, leaving 
probably more than one 5— function contributions to the spectral weight. We then conclude 
that the inclusion of closed loop paths does not suppress the first quasiparticle weight, but 
completely washes out all the quasiparticle excitations at higher energies. Even the few 
peaks that appears in the incoherent part cannot be associated to "string state" resonances 
-as suggested by but are rather similar to the = ones, where localized string 



states cannot exist. We thus confirm the conclusions of Poilblanc. [0 

As far as the energy spectrum E{p) is concerned we obviously found that the lowest 
energy state has a finite quasiparticle weight and has momentum p = (0,0), instead of 
(|, |) as commonly accepted for the t-J model. This is because in the t-J^ model the 

spin fluctuations are neglected, while they play an important role for the energy dispersion 
E{p). 

Finally for the quasiparticle weight as a function of we can apply the same argument 
at the end of the previous section by assuming that Z is basically the overlap of an Sz = 
polaron state of size ^, where ^ may be roughly identified as the correlation length within 
the string picture ^ oc J^^^^- We should get essential singularities like 
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Z oc J-^l^e-^^^ ^'^ for 2C 

-A 

Z oc J7/I^e'^ for d>\. (53) 

As in the string picture, where we could not detect the correct Z Jz behavior for reasonably 
small values of (see Fig.(|^) ), we expect that these kind of singularities are important 
only at a value of ~ 10^^. For larger values of one obtains a crossover to a J^^^ 
behavior surprisingly valid for a quite large range of Jz both in the string picture and in 
the realistic cases shown in Fig.( In the latter picture it is also evident that at some 

small value of Jz the Z factor should vanish much faster than Jf^'^, otherwise we should get 
an unplausible critical value of J^, where the quasiparticle weight vanishes. This at least 
supports the behaviour for Z shown in Eq. (|53D . 



C. Ground state properties for small Jz 

The accurate determination of the ground state energy in the small Jz region is important 

to detect a possible transition between a uniform antiferromagnetic state and a state where 

the hole fully polarize the spins in a finite region of space with size region being somehow 

phase-separated from the remaining antiferromagnetic region. The latter variational state 

has an energy approaching the Nagaoka energy for ^ ^ oo. At finite J^, by optimizing the 

size ^ of the ferromagnetic region, corrections to the asymptotic Nagaoka energy scale like 
2 



Support for a possible transition comes from the infinite dimension limit [|T^. In fact 
the exact evaluation of the Green's function, possible in infinite dimension, implies that the 
lowest energy state contributing to the spectral weight is analytically connected at small Jz 
not to the Nagaoka energy oc —d but to the much higher Brinkman-Rice energy oc 
If the latter state is identified with the lowest energy state in a uniform antiferromagnetic 
phase, there should be a critical value Jc where the 'phase separated' state , discussed at the 
beginning of this section, becomes lower in energy for Jz < Jc- This is essentially the 
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scenario proposed by Emery et al. for the doped t — J model. We will show in the following 
a clear numerical evidence that the above scenario is not confirmed in the t — model for 
a single hole, at least in the 2C case. 

Contrary to the string picture or the infinite dimension limit we see in Fig.(|13|) that the 
asymptotic value for the energy is clearly given by the Nagaoka energy, although the leading 
corrections to the energy seem quite well fitted by the string picture exponent J^^^- In fact 
the diagonal elements of the Hamiltonian in the Lanczos basis describe approximately a 



linear potential (Fig. ([l^) ) as in the string picture. This is obviously important for the small 
Jz correction to the energy. Moreover in Fig.(0) it is shown that the 'phase separated' state 
is well above the estimated energies even for very small J^, leaving a possible transition at 
an unphysically small value of ~ 10^^ for the 2C case. For the 2D case however we do not 
have enough accuracy as shown in Fig.(|D and we cannot exclude a transition at ~ 10~H. 

D. Spin arrangement in the ground state 

Although we have found evidence that the Nagaoka energy is the ground state energy of 
the t-Jz model for — it is not clear what spin background is favoured in this limit. For 
instance we could have that the Nagaoka state is degenerate in the thermodynamic limit 
with an antiferromagnetic state. 

In order to solve this issue, we have calculated the hole-spin-spin correlation function, 
by measuring the one involving the spins in the 2;-direction 

(54) 

and its spin rotation invariant version 

C^iR,) = N< ^lj,\hlSn^Sn^+rMp > ^^^^ 

= < So\SR.SR.+rJSo > 

where is the hole creation operator at origin. These correlation functions measure how 
the spin background is perturbed by the hole. 
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The symmetrized correlation function C^{Ri) has been introduced since, for 0, 
the total spin is a well defined quantum number and the isotropic hole spin-spin correlation 
does not depend on the polarization of the total spin. Thus even in the Sz = sector a 
polaron solution with maximum spin leads to a maximum C{Ri) = |, whereas C^{R) = 
for this polarized state since the contribution of the parallel spins are exactly cancelled by 
the one of the antiparallel spins (all these contributions have the same weight in the polaron 
solution ||29|). In this way, we can unambiguously distinguish the ferromagnetic region with 



C^{Ri) > from the antiferromagnetic one with C^{Ri) < 0. As shown in Fig. |I5|, we have 
found that for large Jz all the spins are antiferromagnetically correlated in the z direction 
since the two previous correlation functions are almost identical. There is a clear evidence 
of an antiferromagnetic correlation which approaches the asymptotic value C{R) —* —1/4 



with a correlation length consistent with the one shown in Fig. (11 



For small Jz, as expected, the finite n corrections are important and tends erroneously to 
enhance the antiferromagnetism (see Fig. (|16|)). Instead by studying the behavior of C{R) 
and Cz{R) as a function of 1/n it is quite clear that the spins are strongly correlated in the 
x — y plane and for n —>■ oo and Jz = 0, the hole-spin-spin correlations seem to approach the 
fully polarized values C{R) = 1/4 and Cz{R) = 0, presumably at any finite distance from 
the hole. 



VI. CONCLUSION 

In conclusion the physical picture of the one hole ground state in the t — Jz model seems 
clear. As we decrease Jz, we approach the Nagaoka state with maximum spin and with 
total Sz = (which is by the way a conserved quantity) due to the proliferation of closed 
loop paths that strongly enhance the ferromagnetic correlations around the hole in the x — y 
plane. Thus there exists a ferromagnetic polaron with a length which diverges for Jz 0, 
but that is continuously connected to the antiferromagnetic length obtained at finite large 
Jz- In this way the polaron state is somehow similar to the phase separated variational 
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state but we have to emphasize that in our approach the locahzed polaron conserves 



the translational symmetry because it is defined after the Gahleo transformation (J^) and 
consequently it is not "phase separated". Moreover the total spin projection on the 2;-axis 



vanishes, as it is conserved locally by the effective Hamiltonian (|1^). The energy of this 
state is much smaller than the "phase separated" variational state and represents a more 
accurate picture of the single hole ground state in the t — model. 

We have presented here a successful attempt to go beyond the retraceable path approxi- 
mation and the string picture for the hole dynamics in an antiferromagnetic spin background. 
A new Lanczos-type of analysis of the Hamiltonian enabled us to get very accurate results for 
the 2C problem and qualitatively similar ones for the 2D case. At = 0, a clear dispersion 
of the main incoherent peak of A{k^ uo) both for the 2C and the 2D case was found. Contrary 
to the prediction of the large spatial dimension we found, at least for the 2C model, that 
the bottom of the incoherent band is dispersionless and coincides with the Nagaoka energy 
e^i i.e. the minimum possible energy by the Nagaoka theorem. This resolve a controversy 
recently proposed by Miiller-Hartmann and Ventura. In fact, based on the Lanczos Spectra 

-A 

Decoding method, we have given an argument implying that Aiuj) oc (tu — eiv)^'^^^^"'"e'""'^''' 
and if for instance A oc {tBR — ^nY, we can easily understand why the infinite dimension 
limit gets no weight for u < ebr. As the dimensionality is increased the tail of the spectral 
weight below the Brinkman-Rice energy vanishes exponentially, though remaining always 
finite up to the Nagaoka energy in any finite d. 

For finite J^, contrary to the string picture, we found only one coherent quasi-particle 
weight and an incoherent broad spectrum at higher energy. A second quasiparticle peak 
may appear in the spectral function but has always a very small weight. The possibility of 
a phase transition as a function of is not compatible with our numerical results for the 
ground state energy, unless for very small coupling constant, and consequently, the Emery's 
argument about the phase separation for small J is found unlikely in the t-Jz model. 

The spin charge decoupling for a single hole is surely not evident for short distance 
propagations. However, it still remains open whether the spin charge decoupling happens 
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asymptotically at large distance, although for the 2C case we have ruled out this possibility 
up to a distance of about 10 lattice spacings. 

We gratefully acknowledge useful discussions with A. Parola, E. Tosatti and D. Poilblanc. 
This work has been supported by CNR under Progetto Finalizzato "Sistemi informatici e 
calcolo parallelo" . 

APPENDIX: QUASIPARTICLE WEIGHT, GREEN'S FUNCTION AND 

CURRENT OPERATORS 

An important quantity to characterize the dynamics of the single hole is the so called 
quasi-particle weight appearing as the residue of a simple pole in the one-hole dynamical 
Green's function. This residue can be alternatively calculated by means of the overlap 



of the ground state \tpp > of one hole with momentum p and the state Cp^a\H >, where \H > 
is the translation invariant ground state without holes: 

Zp=\< H\cUi;p > r = I < H\So > \' (Al) 

where we have explicitly used that n^^o\So >= \So >■ For t > 0, i.e. positive time, the 
Green's function is defined as 

G{p,t) = -1i < H\cl^e-'^^-'^-^''^'cp^^\H > (A2) 

where Eq is the corresponding energy of the state \H >. Here the factor two comes from 
the requirement that G{p,t 0^) = —2 < H\np^cr\H >= —i. The normalized state 
\il)p>= \/2cp^a-\H > is of the form (|^), if we choose 

= V2n^,o\H> . (A3) 

Due to the correspondence of eigenstates between H^^^ and H, we can expand 15*0 > in 
terms of eigenstates of -ffp-^-^ and easily check that the propagation of l^o > with the effective 
Hamiltonian, 15*0 >t= e^^p**^\So >, corresponds exactly to the propagation of ipp with the 
exact t — Jz Hamiltonian and 
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G{p, t) 



(A4) 



Using that \Sq >= \^ni^„\H >, that the commutator H^f^^Ua^o 



vanishes and that G 



does not depend on a, we get, after Fourier transform G(p, cu) = J dt G{p,t) e 





ILdt 



G{p,u) =<H\ 



uj + i6- Hp^^ 



\H > . 



(A5) 



Another important quantity is the current operator, which is useful when we calculate 
the transport properties. On a discrete lattice, it is defined as H: 



(A6) 



where cq is the electron charge. The matrix elements of the current operator between two 
one-hole states with given momentum p define an effective current operator j^^^ acting on 
spin states only: 



<<b.l^>=< s'\fjf\s> 



(A7) 



Analogously to the calculation shown in Section II, the effective current operator can be 
written as 



3 



(A8) 
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FIGURES 

FIG. 1. (a):The ratio of the spectral weight at different sites for the two chain lattice. The 
spectral weight as a function of R was obtained by Fourier transforming A{p,uj) obtained by the 
standard Lanczos Spectra Decoding. Moreover since at fixed n, A(R, u) is exactly zero for R large 
enough, only a finite number of momenta are necessary to implement exactly the mentioned Fourier 
transform (b):The calculated spectral weight for R = O. The solid line is got by summing the 
spectral weight for all momenta verifying the independent relation A(R = 0,lu) = J ^A{p,uj). 
The data points are calculated directly by the Lanczos Spectra Decoding using a trial state with 
a hole localized at the origin, i.e. without using the translation invariance. 

FIG. 2. The application of the effective t Hamiltonian([T^) on the trial state. (a):The Neel 
state with one hole. The hole is located at origin. (b):The state after the action of Tr^. T^^ 
translates the Neel state one lattice space along the direction fi. The hole moves to the nearest 
neighbor. (c):The final state. The spin exchange operator moves the hole back to the origin and 
leaves an over-turned spin defect in the Neel background. 

FIG. 3. The single hole quasiparticle weight of the t-J^ model on the Bethe lattice with z = 3 
and 2 = 4. The inset is an expansion of the small Jz region and the axes have been scaled by a 
factor 1000. The Lanczos matrix was truncated after n = 40000 Lanczos steps, by far enough to 
obtain convergent n = oo results even for very small Jz- 

FIG. 4. The Lanczos density of states on the Bethe lattice (for z = 3 and z = 4), the 2C 
and the 2D lattice. Triangles, squares and circles correspond to the small, medium and large n 
calculation respectively. The continuous lines are the exact results for the Bethe lattice and guides 
to the eye for the 2C and 2D lattices 
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FIG. 5. Plot of the lowest eigenvalues of the 2C and 2D model as a function of 1/n, the inverse 
of the Lanczos-iteration number. For the 2C case , the wavevector k ranges from from (0, 0) 
( bottom) to (vr, 0) (top) with nine equally spaced values. For the 2D case the /c— path in the 
magnetic Brillouin zone is F ^ M ^ X — > F, where F = (0,0) (bottom), M = (7r,0) (top) and 
X = (7r/2,7r/2). The horizontal dashed lines denote the Brinkman-Rice ground state energies. 

FIG. 6. The expected smooth quantity Z[ijo) = (n + 1) x Z at = plotted as a function 
of the energy for different Lanczos iterations. The continuous line connects the n = 26 data. By 
comparing the data at different Lanczos iterations, Z[u) seems to be non vanishing just above the 
Nagaoka energy e^v = — 3t for 2C. The insert is an expansion of the band edge. 

FIG. 7. The linear intercept fj,{n,n — 2) defined in the text ( |48| ) as a function of ^ for 
p = (0,0) and R = (0,0). The dashed lines are linear extrapolations of the last six points. The 
two horizontal arrows indicate the Brinkman-Rice and Miiller-Hartman- Ventura estimates. The 



vertical arrow indicates the largest n analyzed in |20|. 



FIG. 8. Same as in Fig. (0) for the 2D case, the data are taken from [20|. 

FIG. 9. The behavior of Z{uj) at band tail for p = (0,0) and p = ('/r,0). The solid lines are 
the least square fit by expression (^). 

FIG. 10. The spectral function at = 0.3t and k = {tt,0) and Jz = 2t and k = (0,0). Zi and 
Z2 are the quasiparticle weights for the lowest two eigenstates. The small value of Z2 compared to 
Zi indicates that only the lowest weight remains finite for n — > 00, i.e. contributing to the spectral 
weight with a true 5— function. 

FIG. 11. The 2C-wavefunction of the single hole in the Lanczos basis for = 2t, p = (0,0) 
and Jz = 0.3t, p = (7r,0). The solid line denotes the ground state, the dotted one the first excited 
state and the dashed one the second excited state. 
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FIG. 12. Calculated quasiparticle weight Z as a function of {Jz/tf^^- 

FIG. 13. The ground state energy as a function of Jz^^ ■ The data points refers to n = 26 (2C) 
and n = 14 (2D). The continuous lines are a fit = a + bJz + cJz of the data. The dotted line 
is the energy of the "phase separated polaron" described in the text. 

FIG. 14. The diagonal matrix elements of the Hamiltonian in the Lanczos basis. The solid line 
refers to the two-chain lattice and the dashed line to the corresponding Bethe lattice with 2 = 3. 

FIG. 15. The hole spin-spin correlation function C^(i?) and C^{R) with = (1,0). For 
large Jz C{R) is purely antiferromagnetic. However for vanishing (or small) Jz, the ferromagnetic 
component in the x — y plane is important. 

FIG. 16. The extrapolation n ^ oo of C'^{R) and C^{R) for the three lattice sites closest to 
the hole at Jz = 0. The horizontal line is the value of C^{R) for the fully polarized Nagaoka state. 
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TABLES 

TABLE L The non vanishing coefficients Cn{R) of the t-Jz Hamiltonian on the 2C lattice. 
Only those coefficients which are not included in the Miiller-Hartmanni- Ventura's table are shown. 



The notation is similar to the one in r ef. i.e. = \R\^ and n^2 is the number of skeleton paths 



at 


a given distance 


and for a given 


direction (note that there is 


an extra factor two for R ^ O \i 


we 


do not distinguish the two possible directions as 


in ref. Q). 


The data for n = 


26 are accurate 


up 


to ±5. 














2n 


Cn{(f = 0) 


Cn{(f = 2) 


Cn{<f = 4) 


Cn{dP = 10) 


Cn{dP = 16) 


Cn{dP = 26) 


32 


193448 


152314 


108676 


70960 


39081 


17614 


34 


590154 


472488 


340675 


223698 


128823 


58446 


36 


1824844 


1471492 


1068182 


708496 


414196 


196574 


38 


5677040 


4609274 


3385018 


2264848 


1348874 


670086 


40 


17818480 


14539266 


10760828 


7287326 


4397638 


2245908 


42 


56220728 


46154304 


34459409 


23573566 


14437674 


7572526 


44 


178693158 


147425926 


110826815 


76581474 


47561630 


25474418 


46 


570790364 


473551402 


358473393 


249911680 


157303528 


85932926 


48 


1834737522 


1529492974 


1164976270 


819090516 


521903825 


290492088 


50 


5926011194 


4963905566 


3804186739 


2695446152 


1737377480 


983973358 


52 


19240493885 


16187397249 


12476330859 


8905934658 


5800668507 


3338934862 
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2n Cn{d^ = 36) Cn{(f = 50) Cn{d^ = 64) Cn{(f = 82) Cn{d^ = 100) Cn{d^ = 122) C„(d2 = 144) C{p = 0) 



32 5714 1140 984446 



34 


19510 


4294 


534 














3087090 


36 


70917 


18258 


2981 














9727036 


38 


251056 


67790 


12460 


1190 











30898232 
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883914 


258922 


55639 


7634 











98692630 


42 


312G095 


9C)832G 


22232G 


35058 


2GG1 








317324618 
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10934279 


3569044 
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19237 








1025581138 
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472897 


47830 
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464563759 


173622392 


50912034 


11813430 


2117499 


258896 


13277 


35702392358 


52 


1616935261 


626834742 


192138728 


47077158 


9171222 


1332690 


117733 


117646241223 
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TABLE II. 


Fit of the coefficients Cn (R) for R = using 


; the ansatz (pq 


), in the given intervals 


of n given in the rightmost column. 








a 


f^R=o 


Cr=o 


2n 


1.88136 


2.39084 


1.26191 


22-36 


1.88356 


2.42844 


1.38458 


22-38 


1.89515 


2.63760 


2.34793 


22-40 


1.89049 


2.54888 


1.86730 


22-42 


1.89824 


2.70449 


2.81409 
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1.89667 


2.67137 


2.57441 
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1.90428 


2.83985 


4.08250 
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1.90575 


2.87370 


4.48588 
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1.91050 
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6.20875 
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